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We propose a hybrid system with quantum mechanical three-body interactions between photons, 
phonons, and qubit excitations. These interactions take place in a circuit quantum electrodynamical 
architecture with a superconducting microwave resonator coupled to a transmon qubit whose shunt 
capacitance is free to mechanically oscillate. We show that this system design features a three-mode 
polariton-mechanical mode and a nonlinear transmon-mechanical mode interaction in the strong 
coupling regime. Together with the strong resonator-transmon interaction, these properties provide 
intriguing opportunities for manipulations of this hybrid quantum system. We show, in particular, 
the feasibility of cooling the mechanical motion down to its ground state and preparing various 
nonclassical states including mechanical Fock and cat states and hybrid tripartite entangled states. 


Quantum control of macroscopic objects is of great 
fundamental importance [1] and massive mechanical res¬ 
onators strongly interacting with well-controlled quan¬ 
tum systems, e.g. photons and atomic excitations, are 
desired candidates for this purpose W- These can be 
employed for preparing non-classical states in mechani¬ 
cal resonators IMS!, but are also of technological interest, 
e.g. for weak force sensing |10j and transduction of quan¬ 
tum information in quantum networks m- Particularly, 
nonlinear quantum phenomena are very desirable for the 
above purposes as they considerably extend the options 
for manipulation and control of quantum systems. Intro¬ 
ducing an anharmonic part into a setup, can for example 
strengthen its couplings and enriches its physics via the 
nonlinearities [IMS]- Here, we propose a circuit elec¬ 
tromechanical hybrid architecture that combines a nano¬ 
mechanical degree of freedom with both, an intrinsically 
nonlinear component in the form of a superconducting 
qubit and nonlinear interactions between the mechanical 
mode, the qubit excitations, and a harmonic mode of an 
electrical resonator. As a key novelty, this architecture 
features three-body interactions between the mechanical 
and two electrodynamical degrees of freedom, which can 
not be approximated by effective two-body interactions 
due to the involved nonlinearities. 

We explore a circuit quantum electrodynamical sys¬ 
tem consisting of a transmon qubit strongly coupled to 
a superconducting microwave resonator. In addition, 
the resonator-transmon system interacts with a nanome¬ 
chanical oscillator. The advantage of a transmon is its 
robustness against fluctuations of background charges 
achieved by increasing the ratio of Josephson and charg¬ 
ing energies Ej/Eq [U] at the cost of a reduced anhar- 
monicity. Nonetheless its nonlinearity can still be ex¬ 
ploited for controllably producing single photons in a su¬ 
perconducting transmission line resonator via excitation 
exchange m or for controllably producing propagating 


surface acoustic phonons m- Moreover, Josephson junc¬ 
tions integrated into a circuit electromechanical device 
can enhance its optomechanical couplings [El l20H22] . 

Here, we propose to couple the transmon-cavity sys¬ 
tem to a mechanical resonator by replacing one of the 
transmon’s shunt capacitor legs with an oscillating nano¬ 
beam, c.f. Fig. [2 Hereby, we introduce a nonlinear cou¬ 
pling between the qubit and the mechanical resonator 
that can reach the strong coupling regime, i.e., the bare 
coupling rate can exceed the relaxation rates of the sys¬ 
tem. In particular, as the nano-mechanical transmon is 
embedded in a microwave cavity the hybrid system fea¬ 
tures an electromechanical three-body interaction with a 
flux tunable coupling rate. To show some assets of the 
system, we exploit these couplings and the anharmonic- 
ity of the transmon qubit to cool down the system to its 
ground state and prepare it in nonclassical states such as 
mechanical Fock and Schrodinger cat states. 

The model. Our system is composed of a superconduct¬ 
ing coplanar waveguide resonator equivalent to an LcCc 
oscillator, capacitively coupled to a transmon qubit via 
a gate capacitance Cg. The shunt capacitance Cb of the 
transmon qubit depends on the position of a mechanical 
resonator as depicted in Fig. [^b). The Hamiltonian of 
the complete system can be written as {h = 1) [E] 

H = Ho + Hi+Hd, ( 1 ) 

Ho = LOto)a — X{a^)^a^ + Ulrnbrn + cocc'c J- ix(ac^ — a^c). 
Hi = [gt&^a + igtcidc^ - a^c)] {b + b^), 
iJd =fL(ce“^‘ + c1'e-“^*). 

Here, oJt = Ec{^/^—l) with ( = Ej/Ec is the transition 
frequency between the ground and the first excited state 
of the transmon, which is modelled as an anharmonic os¬ 
cillator with annihilation (creation) operator a (a^) and 
Duffing nonlinearity A = Ac/2. The charging energy of 
the qubit is Eq = e^ I2C^ with Cs = Gg J- Cb + Gj (Gj 
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FIG. 1. (color online), (a) Circuit diagram and (b) sketch 
of the hybrid system composed of a transmon qubit, a super¬ 
conducting coplanar waveguide resonator and a mechanical 
oscillator. 


is the capacitance of the Josephson junction). The nano¬ 
mechanical resonator with natural frequency rim is de¬ 
scribed by phononic operators b {b^) and its displacement 
is given hy x = x^.p^(b + &^), where x^pm = 
is its zero-point motion amplitude and to its effective 
mass. The microwave cavity oscillates with frequency 
Wc and is characterized by the bosonic mode operators 
c and o'. An external microwave field drives the cavity 
with amplitude Sl- The rate at which the qubit couples 
to the transmission line is y = 4ifcn.ac(C/2)* with riac 
the rms number of cavity induced Cooper pairs. More¬ 
over, the interactions between the mechanical resonator 
and the other parts of the system are quantified by the 
coupling rates gt = for the transmon-mechanical 

mode interaction and g^c = 45o?^ac(C/2)* for the three- 
body electro-mechanical mode interaction^ where g^ is the 
bare coupling constant [23]. Notably, both couplings are 
enhanced as the ratio Ej /Eq increases, which can be 
tuned in situ and in addition has the beneficial effect of 
enhancing the coherence time of the transmon. Finally, 
we mention that a rotating wave approximation (RWA) 
is applied to get the Hamiltonian H. This is valid for 
C S> 1 and ^ WtjWc, which is compatible with 

the operation regime of our hybrid system. 

In addition to the coherent evolution described by the 
above Hamiltonian, the system is affected by dissipation. 
The energy relaxation rate of the transmon is 7 t = 1 /Ti 
and its total dephasing rate is 1/T| = l/(2ri) -|- l/T^ 
where 7 ^ = is the pure dephasing rate. High qual¬ 
ity superconducting qubits can be fabricated with relax¬ 
ation and dephasing times as high as Ti « 50 /rs and 
T| « 20 /rs [23]. Yet, even higher values, Ti « 70 gs 
and T 2 « 90 ^s, have been realized for 3D cavity setups 
[251126]. In addition, the microwave photons of the cavity 
are subject to loss at a decay rate of Kc and the mechan¬ 
ical resonator is coupled to a thermal bath with the rate 
Tm = where Qm is its mechanical quality factor. 

These dissipation processes are captured by a Liouvillian 
in Lindblad form. Thus, the full dynamics of our system 
is described by the master equation 

P = -i[H, p] + yt'DaP + 14,'Da^aP + (n + l)T^'D~^p 
-t-nT mT’gt P + (2) 

where VoP = opo^ — (6’^dp-|-p6^6)/2 is the dissipator and 


n = [exp(^^) — 1 ] ^ is the thermal phonon number. 

Polariton-mechanical mode interaction. For typical 
configurations, the transmon-cavity interaction is in the 
strong coupling regime (% ^ 7t,«^c)- It is therefore 
convenient to describe the subsystem of cavity photons 
and qubit excitations in terms of dressed state exci¬ 
tations called polaritons, which decouple the interac¬ 
tion zx(acl — die). In terms of these polaritonic modes 
p± = a±d =F iazpC (the explicit forms of a± are given in 
[23]). the Hamiltonian ([^ reads 


Hp = 


X! \^kpl'Pk-Xk{piffk + pnjk^biisx) 


fc=± 




9kpipk + G{p\p-+ p+p^_) (S-fdl), (3b) 

A:—± 


where the polariton resonances are given by uj± = 
a^ujt + Qf^Wc ± 2a+a_x, while the polariton-mechanical 
mode coupling rates are g± = a^gt -\- 2Q!_|_a_gtc and 
G = a+a-gt -\- (a^ — a3)<?tc, and the nonlinearity for 
each polariton is A± = a%X (Note that G > gtc due 
to the contribution from gt). There are thus three-body 
interactions for which the nonlinearity A precludes the 
linearization of the terms p\_p- or p+pl_ . This is in con¬ 
trast to the original three-mode interaction with strength 
gtc or to optomechanical couplings in standard settings 
m- We have included all the inter-polariton interac¬ 
tion terms in H.^ _, see [23]. Note that except for an 

intensity-intensity interaction, all these inter-polariton 
interactions can be neglected in a RWA [23] provided 
UJ+ — uj- = [A^ -|- 4x^]5 A, where A = wt — Wc is the 
detuning between the transmon and cavity frequencies. 
This condition calls for a large photon-qubit coupling 
rate (x S> A) and/or large transmon-cavity detuning 
(A A). To ensure the validity of the above RWA, 
we work in the off-resonance regime with A ^ A. The 
polariton-mechanical interactions in Eq. (3b I provide us 
with a toolbox for the quantum control of the state of 
the mechanical resonator, which is our main interest for 
the proposed architecture in this letter. 

Cooling the mechanical resonator. A first question of 
interest is whether our setup allows for ground state cool¬ 
ing of the mechanical mode as this is a prerequisite for 
many state preparation protocols. We show that this is 
indeed feasible using sideband cooling [25]. In principle, 
both interactions of (3b I are capable of performing the 
task. For the three-mode interaction one would transfer a 
mechanical phonon and a lower polariton excitation into 
a higher polariton excitation, w_ -I- = a;_|_, which sub¬ 

sequently decays. However, a simpler and more efficient 
route is to use the couplings p±p±(b-\-b'^) at large cavity- 
transmon detunings. Here, one polariton is dominantly 
photon-like, while the other describes mainly a transmon 
excitation. In this regime, the photon-like polariton is 
practically decoupled from the mechanical mode, while 
the transmon-like polariton strongly interacts with it at 
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FIG. 2. (color online). Mechanical final phonon numbers ver¬ 
sus drive detuning when only one of the polaritons p+ (solid 
lines) or j5_ (dashed lines) is driven: (a) parameter set^^^^l with 
Q = 150 ± 50 and (b) parameter set^2 with Q — 142 ± 60. 


a coupling rate close to gt- For this kind of interaction, 
the final occupation number of the mechanical mode is 
limited by the total dephasing time T 2 of the qubit and 
the ground state can be reached for > 1 [H]. 

We numerically solve Eq. ([^ with Hamiltonian Q for 
two different sets of parameters. Set#l: Cc = 150 (the 
resonance value, i.e., Wt(Cc) = Wc), m = 1 pg, 50 = 18.2 
kHz, Hm/27r = 10 MHz, Kc/Stt = 10 kHz, 7t/27r = 3 
kHz, and = 8.5 x 10“^. Set#2: Q = 142, m = 3 
Pgj ffo = 20.6 kHz, r2m/27r = 1 MHz, Kc/27r = 50 
kHz, 7t/27r = 5 kHz, and Uac = 1-4 x 10“^. We 
also take 7 ,^ = 27 t and use the common parameters 
Ec/27r = 0.5 GHz, ujJ2t: « 17 GHz, and = 10®. 
The considered mechanical parameters are compatible 
with experimental reported values [5D1132) . In Fig. [^we 
plot the numerical results for the final phonon numbers 
achievable by cooling via either of the polaritons. We find 
that it is possible to cool the mechanical resonator from 
n « 20 phonons (corresponding to an ambient tempera¬ 
ture T ~ 10 mK) to rif « 0.12 for the parameter set#l 
and from n « 100 phonons (environment temperature 
T ~ 5 mK) to ny « 0.03 for the parameter set#2. The 
multiple cooling resonances apparent in Fig. ^b) are a 
signature of the nonlinearity of the coupling [3^ and thus 
provide a measurable witness for the nonlinearity of the 
transmon-viberational mode interaction. Having shown 
the feasibility of ground state cooling, we now describe 
two state preparation protocols enabled by our device. 


Mechanical Fock states. We first describe a protocol for 
preparing the mechanical resonator in Fock states. Our 
strategy here is to first generate individual polaritonic 
excitations and then transfer them to the mechanics via 
the three-mode interaction in Eq. (3b). To this end, side¬ 
band cooling first brings the mechanical resonator close 
to its ground state. Properly shaped microwave pulses at 
suitable frequencies can generate single-qubit rotations 
for the polaritons [M]. Hence, the polariton with higher 
frequency is excited by such a pulse. Then the transmon 
is tuned to the point where ui+ —oj- = and the system 


FIG. 3. (color online), (a) Fidelity of prepared single phonon 
state as a function of the transmon-mechanical coupling rate; 
for ideal initial ground state (solid line) and attainable ground 
state (dots). The point marked by the red circle corresponds 
to the parameter set^l. (b)-(e) Wigner quasi-probability dis¬ 
tribution of the prepared Fock states with grt/^c ~ (b) 5, (c) 
10, and (d) 21. (e) Ideal single phonon state, for comparison. 


evolves for a time ti = 7r/2G, which converts the higher 
energy polariton into a lower energy polariton and a sin¬ 
gle phonon in the mechanical resonator. The generated 
lower energy polariton is finally annihilated by another 
microwave pulse, leaving the system in a single phonon 
Fock state. 

Strong three-body interactions and hence fast excita¬ 
tion transfer could of course be achieved for wt « lOc 
[23] . Yet, in this regime all polariton-mechanical inter¬ 
actions have the same strength =5- = G ~ gt/2^ 
which enables additional undesirable transfer channels 
that hamper the protocol. Hence, we demand for suf¬ 
ficiently large mechanical frequencies to suppress these 
unwanted interactions via a RWA and simultaneously 
ensure u}+ — lu- = Hm. At the same level of accuracy 
the three-mode interaction becomes G{p\'p-h-\-p+f>^_V'). 
Furthermore, the transfer process must be much faster 
than the decoherence rates of the system, therefore, the 
restrictions on the system are max{ 7 t,Fni} gt ^ Hm, 
where F^ « nFm is the effective mechanical decoherence 
rate. The parameter set#l satisfies the above criteria. 

To provide evidence for its feasibility, we numerically 
simulate the protocol, including the initial cooling to the 
ground state, by solving the full master equation (§. 
As very high fidelities for single-qubit gates have already 
been demonstrated, we neglect errors in the polariton ex¬ 
citation and de-excitation steps. The fidelity for a single¬ 
phonon state prepared by sideband cooling followed by 
the above protocol reaches 70% for the parameter set#l, 
see Fig. and could even be enhanced further by in¬ 
creasing the coupling rate and/or starting from better 
ground states, e.g. by employing qubit reset methods 
[35) . To achieve higher number states, the process can 
be repeated until the target state is reached, where the 
interaction times for the swap need to be adjusted to 
Tn = 7r/2y^G with n being the number of mechanical 
phonons which will be obtained at the end of each stage. 
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Note that, although the inter-polariton interactions are 
considerable, they will not play a role in this process as 
long as the total number of polaritons does not exceed 
unity [53]. Once prepared, the state of the mechanical 
resonator can be analyzed by adapting the measurement 
scheme pioneered in |18j . i.e. tuning the transmon to the 
point where — uj- = Om for various interaction times 
and reading out its excitation probability. 

We note that the existence of the cavity mode is crucial 
for the above protocol since the qubit-mechanical inter¬ 
action is not in the form of a state transfer Hamiltonian, 
nor can it be changed to such form by intensely driving 
the (anharmonic) qubit and the original three-mode in¬ 
teraction is not strong enough to be used in this way. 
Yet, the strong coupling of the cavity to the qubit gen¬ 
erates a strong state transfer interaction of three-body 
form between the two polaritons on the one hand and 
the mechanical resonator on the other hand. 

Tripartite hybrid entanglement. We now turn to pro¬ 
pose a protocol for preparing a non-Gaussian tripartite 
entanglement between qubit, cavity, and the mechanical 
resonator in our system. Such states are of interest both 
from fundamental and technical points of view |36j . In 
order to describe the steps for creating them, we go back 
to the original (not dressed) picture of the system and 
consider an effective three-level model for the transmon 
with ground state |0)t and excited states |l)t and |2)t. 

The mechanical resonator is first cooled down to its 
ground state with the transmon and cavity off-resonance. 
By applying a [7r/2]o<-).i pulse, the qubit is prepared 
in a symmetric superposition of ground and first ex¬ 
cited state. Then we let the mechanical resonator inter¬ 
act with the qubit. As the force exerted on the nano¬ 
beam depends on the state of the transmon, such an 
interaction results in a conditional displacement of the 
mechanical system from the origin of the phase space. 
The maximal amount of displacement is 25t/f^m, which 
is achieved when the interaction duration equals half 
the mechanical oscillation period. However, to have 
two distinguishable peaks in the mechanical phase space 
one needs gt ^ Hm, which is not the case in our sys¬ 
tem. This hurdle can be circumvented by applying a se¬ 
quence of Np regularly spaced [7r]o<->.i pulses to the qubit 
with time intervals equal to half of the mechanical pe¬ 
riod. By choosing an odd number of pulses iVp, apart 
from an irrelevant global phase factor, one arrives at 
;^(|0)t|0)c|/3)m + |l)t|0)c|-/3)m) with/3 = {Np + l)gt/n^ 
[371 ED- In the next step, the mechanical resonator 
can be turned into a superposition of odd and even 
cat states, by applying a [7r/2]o<->.i pulse to the trans- 
mon: = |(|0)t|0)c|V’+)m + |l)t|0)c|V’-)m), which 

is already a bipartite qubit-mechanical entangled state. 
Here, A/±= A/± (|/3) ± | — /!)) with the normalization 

factor Af± = [2 ± ^ is an even/odd cat state. 

Now, a single photon in the cavity can be conditionally 



ATp Np 


FIG. 4. (color online), (a) Fidelity of the prepared odd cat 
state versus Np. (b) Probability of finding the transmon in its 
first excited state for the target state ll&tmc) (red triangle), for 
the state resulting from a simulation starting from the ideal 
ground state (blue square), and a simulation starting from an 
attainable ground state (green circle). The parameter set^2 
is used in these simulations. 

produced via a [7r]i<_>2 pulse that flips the qubit from its 
first to its second excited state |l)t —>■ |2)t. Then a flux 
pulse of duration 7r/2-\/2x sets this transmon transition 
in resonance with the cavity. Therefore, the second qubit 
excitation is transferred into the cavity, leading to 

Id/tnpc) = ^(|0)t|0)c|^+)„. + |I)t|l)c|^-)n.). (4) 

The state @ is a hybrid Greenberger-Horne-Zeilinger 
state [33] . Evidently, Idltmc) could also be reduced either 
to even or odd cat states of the mechanical mode by 
performing a post-selection based on the read out of the 
qubit or cavity state. 

To attain macroscopically distinguishable mechanical 
cat states a minimum number of pulses Np is required. 
On the other hand, the realizable Np is limited by the 
decoherence rates of the system. Thus, the parameter 
reg ime allowing for a successful preparation of the state 
(wl) is TT X max{ 7 t,fi„} < ^ < gt. Note also that, 
employing the third level for cloning the qubit excita¬ 
tions as cavity photons is necessary to get the state Q. 
Furthermore, it is desirable to transfer the qubit second- 
excitations into the cavity fast enough to decrease the 
unwanted displacements, so we also demand y ^ flm- 
Since we need to work with the transmon-mechanical 
mode interaction, we chose off-resonance as the appro¬ 
priate working regime, where the three-mode interaction 
is negligible. 

In Fig. I^a) the fidelity of the prepared odd cat state 
is plotted versus the number of applied pulses. Signa¬ 
tures of its generation can be obtained via extracting 
phonon number probabilities with the method exploited 
in Ref. [ID]. The protocol itself can be verified in an 
easier way by measuring the probability of finding the 
qubit in its first excited state or detecting a single pho¬ 
ton in the output of the cavity as a function of the 
number of pulses Np applied to the qubit. Theoreti¬ 
cally, for the ideal state l^ktinc) these probabilities are 
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P(|l)t) = P(|l)c) = Kl-exp{-2(iVp + 1 ) 252 / 0 ^}). 
Fig. I^b) shows P(|l)t) and the probability of finding 
the transmon in its first excited state (or detecting a 
single photon) for the states resulting from numerical 
simulations of the preparation protocol neglecting errors 
in single-qubit gate operations. The coincidence of the 
simulations with the theoretical ideal curve confirms the 
feasibility of the protocol. 
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Supplemental Materials: Quantum State Engineering with Circuit Electromechanical 

Three-Body Interactions 


THE MODEL 


The system we consider is composed of a coplanar waveguide microwave cavity capacitively coupled to a transmon 
qubit whose shunt capacitance depends on the position of a mechanical resonator. The Hamiltonian of the system is 
given by 

H = AEQ{x){n — rig)^ — Ej cos + Ml^Pb + huj'^c'c, (SI) 


where h and ip are the superconducting charge number and phase operators, satisfying the commutation relation 
[ip, fi\ = i and rig is the offset charge of the device which contains both dc and ac contributions. The mechanical 
resonator is characterized by its natural frequency and phonon annihilation (creation) operator b (u'). Its position 
is given by a; = Xzpm(b + 6^) where x^pm = \/h/2mVl^ is the mechanical zero-point motion. And w' is the ‘bare’ 
cavity mode frequency represented by bosonic operator c. 

For large ratios E,j/Ec = C, the anharmonicity of a transmon is not very high and we describe it by expanding the 
cosine term around ip = 0 and keep up to the fourth power in p, 

H = 4:Ec{x){n — Ugf + + hVLmb^b + fiw'c^c. (S2) 

The offset charge rig is induced by an applied dc voltage, intrinsic defects and by the stripline resonators held. The 
latter yields an ac component in rig, so we write Ug = ndc + fiac(c-|-c^), where riac = CgVol2e is the rms number of the 
vacuum induced Cooper pairs. Here, Vq = A/fiw'/2C'c is the root mean square voltage of the resonator’s vacuum held, 
Cc is its capacitance, and Cg the gate capacitance. For a transmon qubit, the n^c can be set equal to zero because 
its inhuence on the energy levels is negligible. Then the Hamiltonian reads 

H = 4Ec{x)n^ + - ^p^ 

-I- 

+4Ecix)nl^ic + c<)^ 

-8ifc(i)^^acn(c + c'^’). (S3) 


We now Taylor expand the charging energy Eq and keep up to the hrst power of x: 

Ec{x) « Ac -f goib + P), 


where we have dehned a bare coupling constant go = -^^Xzpm with do the equilibrium distance between the plates 
of the shunt capacitor with capacitance Cb and Cb = Cg + Cb + Cj the total capacitance. Bosonic annihilation and 
creation operators can be dehned for the transmon such that 


2Ec 


p = 


2Ec 

A, 


(a - a'l'), 
(a -I- d'^). 


The Hamiltonian can now be rewritten in terms of bosonic annihilation and creation operators 

H = huj[d^d — -^{d + d^)^ 

+hn^b^b + Tiuj'^c'c + 4Acna(.(c -I- c^)^ 

-\-ihxid — d^)(c -I- c^) 

-h^{d-d)f{b + P) 

+h^{b + P){c + CY 

+ihgtc{d - a)){b + b'^){c + o') 


(S4) 
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where huj[ = y/8EjEc is the frequency of oscillations between ground state and the first excited state of the transmon 
and hx = 4-Bc^^ac(C/2)^ is the qubit to transmission line coupling rate. Also, the transmon-mechanical hgt = 
and electro-mechanical %c = coupling rates and the three mode interaction rate as hgtc = 4go’T-ac(C/2)^ have 

been defined. The last line of the Hamiltonian includes the coherent drive of the system by a microwave input. Since 
both go and n^c usually acquire very small values, the electromechanical interaction rate gc is negligible in this setup. 

After applying the rotating wave approximation which is valid when Eq 12huj[ (or equivalently when C 3> 1) for 
the Duffing nonlinearity term and X, gt,5c, <7tc ^ w' for the interactions; the Hamiltonian of the system reduces to 
Eq. (1) in the main text with the effective cavity frequency uJc = uj'^ + The parameters listed in Tableshow 

that this approximation is valid indeed. 


POLARITON-MECHANICAL INTERACTION 

Since the transmon-cavity interaction is in the strong coupling regime (x ^ 7t)7c) the subsystem of cavity photons 
and qubit excitations can be described in terms of polaritons. This polaritonic picture can be useful for eliminating 
the major interaction term and diagonalizes some terms of the Hamiltonian. The nonlinear part of the Hamiltonian 
can bring in polariton-polariton interactions which can spoil this simplified picture. However, if the nonlinearity is 
much weaker than the photon-qubit coupling Ec ^ 2/ix and/or the transmon-cavity frequency difference Eq <C 2fi,A, 
we can neglect such inter-polariton interactions in a rotating wave approximation. For now, we will keep all the terms 
for completeness, however, in some situations we will focus on this regime. We define the following polaritonic modes 


p_i_ = Qf+a — ia-C, 

P- = a-d + ia+c, 


where and Q!_ are real positive numbers satisfying -I- = 1 and given by 


a± 



A 

VA2 -h dx^ 


1 

2 


Finally, we arrive at the following polariton-mechanical Hamiltonian: 


(S5) 


— = -I- - x+p\_p\_p+p+ + UJ-P^_P- - X-P^_P^_P-P- -^ 

+ i9+P+P+ + g-P-P- )ib + P)+ G{p'Ip- + p+pl )ib + b^) 

-a_£L(p+e“^* -f -b a+£:L(P-e“^‘ (S 6 ) 


Here, the following coupling constants have been introduced: Conventional polaritomechanical coupling rates 5 + = 
a\gt + + ^.a+a-ptc and g- = a^gt + a+ffc — 2a+Q;_5tc, and the three-mode interaction G = a+a-{gt — Pc) + 

— a^)gtc- The nonlinearities of the polaritons are A+ = Eccy.\/2h and A_ = Eca^/2h. The inter-polariton 
interaction terms are included in _ 


H+- 


Ec 


ctXp+P^+P+P+ + a'!_pLpLp-P- + 4:aXct^p\-P+P'-P- 


+«+«?. (+P+P+P-P- +p+p+pXpX) 
+2aXa-{pl_P+P+P- +p+p+p+pX) 
+2a+a^_{p+pLpXp- + p\-pXp-P-) 


(S7) 


To have an illustration of the coupling rates, we plot the polariton-mechanical coupling rates in Fig. SI The 
couplings Pc and ptc are not shown in the plot, because the corresponding lines almost coincide with the horizontal 
axis. One notes that on-resonance all the three polariton-mechanical interactions are the same and equal to half a 
qubit-mechanics interaction. We will use this property of the polaritonic interactions to prepare mechanical Fock 
states. Also, for large cavity-transmon detunings the coupling rate of the three-mode interaction becomes negligible 
and can be disregarded in this regime. 
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FIG. SI. Rates of the coupling between polaritons and the mechanical resonator normalized to the transmon-mechanical 
coupling rate gt- 


COOLING THE MECHANICAL RESONATOR 


The interactions between the polaritons and the mechanical resonator make it possible to cool the mechanical 
resonator toward its ground state. To this end we work in the regime where the inter-polariton interactions can be 
neglected via a rotating wave approximation. This approximation will hold for a;+ — a;_ = [A^ +4%^] 2 ^ Ec/2h. For 
the parameters listed in Table|l]this approximation does hold only if A ^ Ec/2H. Working in the off-resonance regime 
means that only the qubit-like polariton is effectively interacting with the mechanical resonator. This is because gc 
and gtc are much smaller than gt, see Fig. |S1[ In such situations, one only drives the interacting polariton mode. 
Therefore, when the cavity and the transmon are off resonance one arrives at the following Hamiltonian for the system 
in the frame rotating at the frequency of the input drive wl: 


— = - 5+p^^p+ - X+pI_p^^p+p+ - 5-pip- - X-p^_pLp-P- 

+{9+P+P+ + 9-P-P-){b + S'l’) + Ao pl_p+pip- 

-a-£L{p++P+) + a+£L{p-+P-)- (S8) 


Here, HAq = 2Eca^a- is the polaritonic intensity-intensity interaction and = wp — uj± is the detuning of the 
polariton. Out of resonance, only one of the polaritons is ‘truly’ fed by the input drive, because optimal sideband 
cooling happens for and this leads to Jzp making one of the polaritons an idler, see Fig. [S2 



FIG. S2. (color online). Avoided crossing and the blue and red sidebands of the polaritonic modes (green and violet bold lines) 
in the plane of Ej / Eq and the normalized detuning of the input drive with respect to the transmon frequency. 
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TABLE I. Parameters of the system. 


Quantity 

set^l 

set #2 

c 

15(0t5O 

142‘*±60 

m 

1 Pg 

3pg 

go 

18.2 kHz 

20.6 kHz 

gtiCc) 

315 kHz 

350 kHz 

x(Cc) 

31.4 MHz 

510 MHz 

Hm/2-7r 

10 MHz 

1 MHz 

n 

to 

10 kHz 

50 kHz 

7t/27r 

3 kHz 

5 kHz 

SEcriac/hu^ 

2.0 X lO""^ 

3.4 X 10“® 


resonance value, i.e., tt)t(C) = 

prefactor of the qubit-cavity coupling rate whose fundamental upper bound is ~ 0.1. 


STATE PREPARATION PROTOCOLS 


Mechanical Fock states 


Here we employ the three-mode interaction to convey the excitations to the mechanical resonator. According to 
Fig. |S1| the regime in which such interactions are significant is when the cavity and the transmon are in resonance 
(wt « u>c)- Moreover, one needs to have |a;+ — a;_| « to make the state transfer. Specifically, we consider the case 
of w_|_ — = Oni which for around resonance means that x ~ ^m/2 must hold, see Fig. |S3[a). 

The following recipe can be used to prepare Fock states in the mechanical resonator: (i) Prepare the system in its 
ground state: |0-|_, 0_, Om)- (ii) Flip up the ‘-h’ polariton by applying a [7r]_|_ pulse to get |l+,0_,0m). (iii) Let the 
system evolve for a duration of ti = 7r/2G'. Assuming fast enough single qubit gate operations this leads to the state 
|0+, 1_, Im)- Note that basically the inter-polariton interactions cannot affect this process because from Eq. (S7) it 
can be easily proven that the sates |0+, 1_) and |l+,0_) do not change under such interactions. However, to produce 
higher Fock states one needs to get rid of the excitation in the ’ polariton mode. Hence, the fourth step of the 
protocol must be: (iv) By applying a [7r]_ pulse the ’ polariton is flipped down |0+,0_,li„)- To reach the target 
state: (v) Repeat steps (ii) to (iv) by replacing ri with = 'njly/nG^ where n is the number of the rounds the 
protocol has been repeated (or equivalently, the number of phonons in the next Fock state). To detect or verify this 
state, one could reverse the steps and measure state of the qubit at the end. 


Tripartite hybrid entanglement 


A tripartite entanglement between components can be prepared in the considered system. We work with the 
original (not the dressed) picture of the system to describe the steps for creating such states. We also consider an 
effective three-level model for the qubit with |0)t its ground state, |l)t and |2)t its first and second excited states. 
Furthermore, because the protocol is mostly executed in the off-resonance regime, the qubit-mechanical interaction 
is the only important interaction of the mechanical resonator while the optomechanical and three-mode interactions 
are negligible. We assume that the cavity is fabricated such that the frequency of its fundamental mode matches LOyi 
of the transmon, that is, when on resonance the cavity shares excitations with the first and second excited states of 


the qubit, see Fig. S3'c). 


The steps for preparing the state are the following: (i) When the qubit and the cavity are off resonance, the 
system is prepared in its ground state |0t, Oc, Om)- (ii) A [7r/2] pulse is applied to the qubit to make the superposition 
state ^ (|0t, Oc, Oni) -b |lt,0c,0m))- This superposition state in the qubit results in a state dependent force on the 
mechanical resonator such that when the qubit is in its ground state the mechanical resonator will experience no 
force, and therefore no displacement and relevant phase shift, while when the qubit is excited a force will be applied 
to the mechanics, (iii) Let the mechanics and the transmon interact for a time interval t'. Omitting a global irrelevant 
phase factor one gets ^[|0t,0c,0m) + ^ |lt, Oc, Q:(t'))] i where the mechanical part is in a coherent state with 

— 1) and the phase difference is — sinHmt')- The relation for a(t') already shows 
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FIG. S3, (a) Resonance frequencies of the system with the dressed resonance frequencies appropriate for the preparation of 

mechanical Fock states, (b) Pulse sequence applied to the qubit for maximally displacing the mechanical resonator with 7 flips 
(solid line). The oscillation of the mechanical resonator is displayed for comparison (dashed line), (c) Energy level diagram of 
the transmon when on-resonance with the cavity, suited for generating the tripartite entangled state. 


that since gt <C fim it is impossible to have a superposition of macroscopically distinguishable mechanical states after 
this single preparation step. However, it is possible to increase the ‘distance’ between the two mechanical states by 
applying a sequence of [tt] pulses at time intervals equal to one half mechanical period (<' = Tr/Hm) [Fig. [S^b)]. In 
this way, by applying Np of such pulses we arrive at ^ [jOt, Oc, —/3o) + Oc, /3i)] where 

/3o-/3i = 2(7Vp + l)^ (S9) 

is the distance between the coherent states and 9 is the relative phase, which is zero for odd number of pulses and 
2 

equals to ^ = 7r(a;t — ^)/^m for even A^p. In order to have a symmetric displacement around origin of the phase 
space, it is reasonable to choose an odd number of such pulses Np which results in ^ [jOt, Oc, (3) + jit, Oc, —/?)] with 
P = (Np -b l)gt/flm- (iv) By applying another [7r/2] pulse to the transmon, we push the state into a superposition 
of even and odd cat mechanical states ^(jOt,Oc, V’-i-) + |lt)Oc,"0-)), where Af±tjj± = A/±(|/3) ± | — /?)) is an even or 

odd cat state with the normalisation factor Af± = [2 zb ^. (v) A [7r]i^2 pulse flips the qubit |l)t —)■ |2)t. 

(vi) A flux pulse of duration 7r/2-\/2x brings the relevant transition of the transmon into resonance with the cavity 
(wi 2 = Wc) and transfers the qubit second excitations to the cavity. Finally we thus arrive at the following state: 

i(|0t,0c,V’+) + |lt,lc,V’-))- (SIO) 

This is a hybrid Greenberger-Horne-Zeilinger state. 





































